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Abstract. Collisions of actinide nuclei form, during very short times of few zs (1CT 2 ' s), the heaviest ensembles 
of interacting nucleons available on Earth. Such collisions are used to produce super-strong electric fields by the 
huge number of interacting protons to test spontaneous positron-electron pair emission (vacuum decay) predicted 
by the quantum electrodynamics (QED) theory. Multi-nucleon transfer in actinide collisions could also be used 
as an alternative way to fusion in order to produce neutron-rich heavy and superheavy elements thanks to inverse 
quasifission mechanisms. Actinide collisions are studied in a dynamical quantum microscopic approach. The 
three-dimensional time-dependent Hartree-Fock (TDHF) code tdhf3d is used with a full Skyrme energy density 
functional to investigate the time evolution of expectation values of one-body operators, such as fragment position 
and particle number. This code is also used to compute the dispersion of the particle numbers (e.g., widths of 
fragment mass and charge distributions) from TDHF transfer probabilities, on the one hand, and using the Balian- 
Veneroni variational principle, on the other hand. A first application to test QED is discussed. Collision times in 
238 U+ 238 U are computed to determine the optimum energy for the observation of the vacuum decay. It is shown 
that the initial orientation strongly affects the collision times and reaction mechanism. The highest collision times 
predicted by TDHF in this reaction are of the order of ~ 4 zs at a center of mass energy of 1200 MeV. According 
to modern calculations based on the Dirac equation, the collision times at E cm > 1 GeV are sufficient to allow 
spontaneous electron-positron pair emission from QED vacuum decay, in case of bare uranium ion collision. 
A second application of actinide collisions to produce neutron-rich transfermiums is discussed. A new inverse 
quasifission mechanism associated to a specific orientation of the nuclei is proposed to produce transfermium 
nuclei (Z > 100) in the collision of prolate deformed actinides such as 232 Th+ 250 Cf. The collision of the tip of 
one nucleus with the side of the other results in a nucleon flux toward the latter. The probability distributions for 
transfermium production in such a collision are computed. The produced nuclei are more neutron-rich than those 
formed in fusion reactions, thus, leading to more stable isotopes closer to the predicted superheavy island of 
stability. In addition to mass and charge dispersion, the Balian-Veneroni variational principle is used to compute 
correlations between Z and N distributions, which are zero in standard TDHF calculations. 



1 Introduction 

Actinide collisions are important tools to test our under- 
standing of the nuclear many-body problem. They form 
nuclear systems in extreme conditions of mass and isospins. 
The prediction of the outcome of such collisions is a great 
challenge for nuclear theorists. 

In particular, the question of "How long can two ac- 
tinides stick together" is of wide interests. The quantum- 
electrodynamic (QED) theory predicts that spontaneous pairs 
of e + + e~ may be emitted due to the strong electric fields 
produced by the protons [1 2 3 1 . This process is also known 
as "QED vacuum decay". It occurs when an empty elec- 
tron state dives into the Dirac sea. QED predicts that such 
a hole state is unstable and decays by e + + e~ pair pro- 
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duction. The life-time for such process is, however, longer 
than the collision-time between actinides. Then, the lat- 
ter has to be optimized to allow for an experimental ob- 
servation of vacuum decay. Recent calculations based on 
the time-dependent Dirac equation [3 1 show that two bare 
238 U need to stick together during at least 2 zs to allow 
for an observation of spontaneous positron emission. Al- 
though no pocket exists in the nucleus-nucleus potential of 
this system [4 5 6 7|, nuclear attraction reduces Coulomb 
repulsion and dissipation mechanisms such as evolution of 
nuclear shapes may delay the separation of the system [8|. 
In a recent experiment, delay times in this reaction was 
searched analyzing kinetic energy loss and mass transfer [9] 

Another application of actinide collisions is to form 
neutron-rich heavy and superheavy nuclei by multi-nucleon 
transfer H 1 01 1 1181 . Such reactions could be used to explore 
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the "blank spot" between decay chains of nuclei formed by 
"hot" and "cold fusion" around Z = 105 and N = 160. 

Theoretically, the complexity of reaction mechanisms 
and the high number of degrees of freedom to be included 
motivate the use of microscopic approaches. Early dynami- 
cal microscopic calculations of 238 U+ 238 U were performed 
with spatial symmetries and simplified effective interac- 
tions 0121131 . Recently, this system has been studied within 
the Quantum Molecular Dynamics (QMD) model 061141 in 
which nucleon wave functions are constrained to be Gaus- 
sian wave packets and with the time-dependent Hartree- 
Fock approach which overcomes this limitation 1 15 1. 



2 Formalism 

In general, the full quantum many-body problem cannot 
be solved exactly and, in most realistic cases, approxima- 
tions have to be made. In general, variational principles 
are useful to build approximation schemes by reducing the 
variational space. 

2.1 The Balian-Veneroni variational principle 

The Balian and Veneroni (B V) variational principle is based 
on the action lfl6l 

S BV = Tr[£>(t i)B(fi)] - f t Tr(fi^ - iD[H,B]^, (1) 

where B and D are the time-dependent trial observable and 
density matrix of the trial state, respectively, and Tr de- 
notes a trace in the Fock space. Both the state and the 
observable are allowed to vary between to and t\, corre- 
sponding to a mixture of the Schrodinger and Heisenberg 
pictures. They are constrained to obey the mixed bound- 
ary conditions Z)( f o) = A) and B{t\) = X, where Dq is 
the density matrix of the initial state of the system and X 
is the operator we want to evaluate at time t\ . Without re- 
striction on the variational spaces, the variational princi- 
ple 6S bv = 0, with the above conditions, is fully equiva- 
lent to the Schrodinger equation if the initial state is pure 
O = \¥ )(¥ \). 



2.2 Mean-field approximation 

In most practical applications, mean-field models are con- 
sidered in a first approximation. In mean-field theories, the 
interaction between the particles is replaced by a one-body 
mean-field potential generated by all the particles. It is then 
assumed that each particle evolves independently in this 
potential. 

For instance, N independent fermions may be described 
by a Slater determinant \(p) = Y[f=i®l\~)y where a\ cre- 
ates a particle in the state \<pi) when it is applied to the 
particle vacuum |— ). In such a Slater determinant, all the 
information is contained in the one-body density-matrix 
p = Xi=i IViKfil The BV variational principle is usually 



applied at the mean-field level where the variational space 
of D is restricted to independent particle states, i.e., with 
D = |0)<0|. 

2.3 Expectation values of one-body operators 

In addition to the mean-field approximation, the variational 
space for B is usually constrained to belong to the same 
class of operators as the observable of interest. For in- 
stance, if one wants to predict expectation values of one- 
body observables X = zZi=\ Qx(f), then it is natural to re- 
strict the variational space for B to one-body operators. In 
this case, one recovers the TDHF equation 0171181 

ih% = [h[plp], (2) 
ot 

where h[p] is the Hartree-Fock (HF) single-particle Hamil- 
tonian with matrix elements h a p — , H is the full 

Hamiltonian, andp ff/3 = {ip a \p\<f>p> = (0|a^2„|0). 

According to this variational approach, TDHF is an op- 
timized mean-field theory to describe expectation values of 
one-body observables. However, TDHF may fail to repro- 
duce their fluctuations a xx = V<£ 2 > - (X) 2 Jl9l20l . 

2.4 Fluctuations of one-body operators 

The B V variational principle can also be used with the vari- 
ational space B e {e 7af "} to determine an optimum mean- 
field prediction for correlations ctxy and fluctuations <Txx 
of one -body operators 0211221 . with 

(Txy = ^(XY) - (X)(Y). (3) 
In case of independent particle states, this leads to 

oi F (f!) = lim -i-tr([jo(fo) - px(to, e)] [p(t ) - p Y (t , e)]) , 
e->o ze 

(4) 

where tr denotes a trace in the single-particle space. The 
one-body density matrices px(t, e) obey the TDHF equa- 
tion (0 with the boundary condition 

Px(h,e) = e^p(t{)e-^\ (5) 

while p(t) is the solution of Eq. (f2| with the initial con- 
dition Papih)) - Tva'^a a Do = (0o|a^o#o)- The optimum 
mean-field prediction of ctxy in Eq. (0]i differs from the 
"standard" TDHF expression which is evaluated from Eq. (0 
using p(ti). 

Eq. (HJi has been solved numerically in the past with 
simple effective interactions and geometry restrictions [23 241251 . 
Modern three-dimensional TDHF codes with full Skyrme 
functionals 02612712 8 29 1 can now be used for realistic ap- 
plications of the BV variational principle 03OI311 . 

In this work, the fluctuations cr^y, <x zz , and <taa, are 
computed in fragments resulting from actinide collisions. 
The correlations cr NZ , which are strictly zero in standard 
TDHF calculations, are also determined with this approach. 
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Fig. 1. Isodensities at half the saturation density, i.e., po/2 = 0.08 far , in 238 u+ 238 U central collisions at a center of mass energy 
£„„. = 1200 MeV. Evolutions associated to the four initial configurations xx, yx, yy, and zy are plotted in columns (time runs from top 
to bottom). Consecutive snapshots are separated by 1.125 zs. 



3 numerical details 

The use of a three-dimensional TDHF code with a full 
Skyrme energy-density-functional (EDF), modeling nuclear 
interactions between nucleons and including spin-orbit in- 
teraction [32 33 1, allows for a realistic prediction of these 
quantities. The TDHF equation (|2]) is solved iteratively in 
time, with a time step At — 1.5 x 10~ 24 s. The single- 
particle wave-functions are evolved on a Cartesian grid of 
96 x 32 x 32/2 points with a plane of symmetry (the col- 
lision plane) and a mesh-size Ax = 0.8 fm. The initial dis- 
tance between collision partners is 22.4 fm. The tdhf3d 
code is used with the SLy4a! parameterization 11261 of the 
Skyrme EDF, which is the only phenomenological ingre- 
dient, as it has been adjusted on nuclear structure proper- 
ties [33 1 . Ref. 1 34 1 gives more details of the TDHF calcu- 
lations. The numerical details for the evaluation of Eq. © 
can be found in |f3T|. 



4 Collision time in 23S U+ 238 U 

TDHF calculations have been performed to investigate the 
collision time in 238 U+ 23S U Q3. The 238 U nucleus exhibits 
a prolate deformation with a symmetry axis in its ground 
state. The effect of this deformation on collision is inves- 
tigated in four configurations (xx, yx, yy and yz) associ- 
ated to different initial orientations. The letters x, y and z 
denote the orientation of the symmetry axis of the nuclei 
which collide along the x axis [see Fig. [Q. We focus on 
central collisions as they lead to the most dissipative reac- 
tions with the longest collision times. 

Here, the collision time T co n is defined as the time dur- 
ing which the neck density exceeds po/10. It is shown in 




250 1500 
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Fig. 2. Collision times for each orientation as function of center 
of mass energy. The shaded area indicates the limit of 2 zs above 
which vacuum decay is expected to be observable in central col- 
lisions. 



Figure [2] as function of the center of mass energy E cm . At 
E c .m. ^ 900 MeV, three distinct behaviors between the xx, 
yx and yy/yz configurations are seen. In particular, the last 
need more energy to get into contact as the energy thresh- 
old above which nuclear interaction plays a significant role 
is higher for such compact configurations. 

At all energies, the yx, yy and yz orientations exhibit 
roughly similar behaviors, i.e., a rise and fall of T co u with 
a maximum of 3 - 4 x 10~ 21 s at E CM . ~ 1200 MeV. Dy- 
namical evolution of nuclear shapes in these three config- 
urations and a strong transfer in the yx one (see next sec- 
tion) are responsible for these rather long collision times 
as compared to scattering with frozen shapes of the reac- 
tants [8|. The xx configuration, however, behaves differ- 
ently. For 700 < E c . m , < 1300 MeV, T co u exhibits a plateau 
which does not exceed 2x 10~ 21 s. This overall reduction of 
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T co u in the xx case is attributed to the strong overlap of the 
tips, producing a density in the neck higher than po [ 1 5 1 . 
The fact that nuclear matter is difficult to compress trans- 
lates into a strong repulsive force between the fragments 
which decreases their contact time. This phenomenon is 
also responsible for the fall of collision times in the other 
configurations, though higher energies are needed to strongly 
overlap. 

The calculations of Ref. 1 3 1 show that the observation 
of spontaneous emission of e + + e~ needs a contact time of 
at least 2 zs between the bare uranium nuclei. The TDHF 
calculations of Fig. [2] predict that such contact times are 
reached in central collisions for energies E cm > 1 GeV. 
This lower energy limit should be taken into account in 
future experimental programs dedicated to the search of 
QED vacuum decay. Finally, it is worth mentioning that 
other approaches lead to comparable collision times in ac- 
tinide collisions B6I35I . 

5 Formation of neutron-rich transfermium 
nuclei 

5.1 238 U+ 238 U reaction 

We now analyze the proton and neutron numbers of the 
fragments produced in exit channels of actinide collisions. 
Strictly speaking, these fragments should be considered as 
primary fragments as they might decay by statistical fis- 
sion. This decay is not studied here as it occurs on a much 
longer time scale than the collision itself. The importance 
of initial orientation on reaction mechanism is clearly seen 
in Fig. Q] for the 238 u+ 238 U reaction. For symmetry rea- 
sons, the xx, yy, and yz configurations give two symmetric 
distributions of fragments, although nucleon transfer is still 
possible thanks to particle number fluctuations. Nucleon 
transfer is expected to be stronger in the yx configuration 
because, in addition to fluctuations, no spatial symmetry 
prevents from an average flux of nucleons. The yx config- 
uration is then expected to favor the formation of nuclei 
heavier than 238 U. 



5.2 232 Th+ 250 Cf reaction 

Similar calculations have been performed on the system 
232 Th+ 250 Cf ll36l . The same effect is observed, i.e., an im- 
portant multi-nucleon transfer in the xy and yx configu- 
rations. The xy configuration where the 250 Cf nucleus re- 
ceives nucleons (its deformation axis is perpendicular to 
the collision axis while the one of 232 Th is parallel to the 
collision axis) corresponds to an "inverse quasifission" mech- 
anism due to a specific orientation of the collision partners. 
Indeed, contrary to standard quasifission, the exit channel 
is more mass asymmetric than the entrance channel. Note 
that inverse quasifission may also occur due to shell effects 
in the exit channel (8) . 

The effect is illustrated in Fig. [3] where the distribu- 
tion of heavy fragments is shown at E cm = 916 MeV. This 
distribution is computed after a TDHF calculation, using 



a particle number projection technique [37 1. The center of 
the distribution is located around 265 Lr, i.e., in the neutron- 
rich side of the known Lawrencium isotopes. Note that, at 
the end of the TDHF calculation, the decay of the frag- 
ments by neutron emission is only partial Il36l . 

The width of such a distribution is known to be un- 
derestimated at the TDHF level II 191201 . In addition, we 
see in Fig. [3] that the probability distributions for N and 
Z are uncorrelated at the TDHF level. This is not a fea- 
ture of the TDHF formalism itself, but rather a limitation 
due to the fact that, for practical applications, one assumes 
\<P) - I0p) x l0n)> where \<p p ) and \<p n ) are Slater determinants 
of the proton and neutron single-particle wave-functions, 
respectively ||3T1 . If this constraint is released, as in [38 1, 
then non-zero correlations could be obtained at the TDHF 
level. 

The Balian-Veneroni variational principle can be used, 
at the mean-field level, to optimize both widths of proton 
and neutron distributions as well as their correlations. Re- 
alistic calculations have been performed recently to study 
deep-inelastic collisions [31 1. Similar calculations have been 
done to investigate the inverse quasifission mechanism dis- 
cussed above. Preliminary results are shown in Fig.|4]where 
the heavy-fragment distribution is shown for the same xy 
orientation of 232 Th+ 250 Cf as in Fig. [3] As expected, much 
larger widths than in the TDHF case are observed in one 
hand, and, in the other hand, strong correlations between 
the proton and neutron number distributions are observed, 
which can be seen by the fact that the fragments are pro- 
duced along the valley of stability. 



6 Conclusions 

To conclude, this fully microscopic quantum investigation 
of actinide collisions exhibits a rich phenomenology strongly 
influenced by the shape of the nuclei. Two main conclu- 
sions can be drawn, (z) The giant system formed in bare 
uranium-uranium central collisions is expected to survive 
enough time with an energy £„„ > 1000 MeV, for the 
spontaneous positron emission to occur, (if) The primary 
heavy-fragments produced by multinucleon transfer are more 
neutron-rich than in fusion-evaporationreactions. The width 
of these distributions, computed with the Balian-Veneroni 
prescription, are much larger than with TDHF. Associated 
cross-sections need to be determined to estimate the exper- 
imental possibility of neutron-rich transfermium and SHE 
productions. Extension of the formalism need to be inves- 
tigated. For instance, the role of pairing could be studied 
with Time-Dependent Hartree-Fock-Bogoliubov (or BCS) 
codes [39 40 41 1 . Stochastic-mean-field methods might also 
be applied to investigate the role of initial beyond-mean- 
field correlations on fluctuations H2l . 
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Fig. 3. TDHF heavy fragment probability distribution (linear scale) for 232 Th+ 250 Cf central collision in the xy configuration. 
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Fig. 4. Same as Fig. [3] with the BV prescription. 
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